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We examine the long-term correlations and multifractal properties of daily satellite retrievals of 
Arctic sea ice albedo and extent, for periods of ^ 23 years and 32 years respectively. The approach 
harnesses a recent development called Multifractal Temporally Weighted Detrended Fluctuation 
Analysis (MF-TWDFA), which exploits the intuition that points closer in time are more likely to be 
related than distant points. In both data sets we extract multiple crossover times, as characterized 
by generalized Hurst exponents, ranging from synoptic to decadal. The method goes beyond 
treatments that assume a single decay scale process, such as a first-order autoregression, which 
(— I ■ cannot be justifiably fit to these observations. Importantly, the strength of the seasonal cycle 

pi J "masks" long term correlations on time scales beyond seasonal. When removing the seasonal 

cycle from the original record, the ice extent data exhibits white noise behavior from seasonal to 
bi-seasonal time scales, whereas the clear fingerprints of the short (weather) and long ('^ 7 and 9 
year) time scales remain, the latter associated with the recent decay in the ice cover. Therefore, 
long term persistence is reentrant beyond the seasonal scale and it is not possible to distinguish 
* ^ ■ whether a given ice extent minimum/maximum will be followed by a minimum/maximum that is 

larger or smaller in magnitude. 
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1. Introduction 

Earth's polar oceans are viewed as a bellwether of climate change because their surfaces 
are covered by a thin (several meters) mosaic of high albedo sea ice floes that modulate the 
atmosphere/ocean heat flux. Whence, as opposed to the massive meteoric ice sheets that 
are several kilometers thick, sea ice is considered to be a more sensitive component of the 
cryosphere to perturbations and feedbacks, particular ly the ice-albedo feedback which has 
driven large scale climate events over Earth history feal tzmanl. l2002h . Indeed, the retreat 
of Arctic sea ice coverage during recent decades (Fig. [T|) has captured substantial interest 
dKwok and Untersteineil . l20Tll : ISerrez3 . l201lh . An essential question concerns the nature 



of the decay in ice coverage; is it a trend associated with the influence of greenhouse 
forcing, or is it a fluctuation in a quantitative record that is short (~ 30 years) relative to 
the dynamics of the cryosphere on climatic epochs (> 10^ yrs)? Our goal is to extract 
the longest intrinsic time scale from the data to provide a framework for geophysical 
modeling. 
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Both past climate data and basic physical arguments indicate that a sufficiently 
large increase in greenhouse gas concentration will drive the decay of the ice cover, 
whereas the state of the art glob al climate models under-project the observed recent decay 
dKwok and Untersteinej . l20Tll) . The stabilizing respo nse to deleterious perturbations of 
the ice cover was examined by iTietsche et al.l (1201 ih . who numerically prescribed ice- 
free summer states at various times during the projection of 2P' century climates and 
found that ice extent typically recovered within several years. Suc h rapid response times 
can b e captured within the framework of relatively simple theory ('Moon and WettlaufeiJ, 



20111) . but independent of model complexity, both internal and external forcings and their 



intrinsic time scales manifest themselves in large scale observations of the geophysical 
state of the system. Due to the fact that we cannot a-priori exclude the observed decline 
in the ice cover as being an intrinsic decadal oscillation or nonstationary influence in 
the climate system, we use the finest temporal resolution in the observed record (days) 
to examine the action of multiple scales, from that of the weather and beyond. The 
fingerprints of the noisy dynamics of the system on time scales longer than the seasonal 
record may reside in that record itself, and our goal here is to extract these to provide 
a framework for geophysical modeling of the underlying processes that created them. It 
is hoped that the approach will provide theoretical constraints as well as guidelines for 
comparative data analysis studies. 

Most observational studies of the satellite recor ds of ice coverage extrapolate in time 
the annual or monthly means (e.g., see Fig. 1 of ISerrez^ . 1201 li ). Whilst the observed 
declines over this troika of decades (particularly the last decade) are striking, our goal 
here is to begin a systematic effort to examine the many processes that influence 
the ice cover viz., the intrinsic time scales reflected in geophysical data. The general 
methodological framework we employ allows one to examine time series data in a 
manner that can distinguish between long term congelations and trends. We examine 
whether there exists a multiplicity of persistent scales in the data that can provide a 
basis for examining cause and effect in the geophysical scale observables of the system. 
Therefore, we view several types of large-scale sea ice data as a multifractal system 
in which a sp ectrum of scaling-the singularity spectr um-characterizes the behavior of 
nearby points dStanley and Meakinl.ll988HFeded.ll988h . The basic approach of relevance 
is the multi-fractal generalization of Detrended Fl uctuation Analysis (DFA) a dvanced 
and widely applied by Kantelhardt and colleagues dKantelhardt et al.L |2002| . 12006 ) aptly 
called Multifractal Detrended Fluctuation Analysis (MF-DFA). In the last decade this 
approach has been developed in many directions, f rom s tudying extreme events with 
nonlinear long term memory (Bogachev andB unda. 201 ih. to examini ng the influence 



of additive noise on long term correlations (e.g.. iLudescher et"an. 1201 lb . Here we use a 
new extension of this methodology called Multifractal Temporally Weighted Detrended 
Fluctuation Analysis (MF-TWDFA), which exploits the intuition that in any time series 
points closer in time are more likely to be related than distant points and can provide a 
ra ther more cl ear signature of long time scales in the fluctuation function and its moments 
(iZhou and Leu ng. 2010). In the next section we summarize MF-TWDFA and describe 
its connection to the conventional approaches of analyzing autocorrelation functions and 
power spectra. The geophysical data are described in Section [3] after which we present the 
principal results from MF-TWDFA in Section |4]bef ore concluding. 
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Figure 1 . Equivalent Ice Extent (EIE) during the satellite era, with the seasonal cycle removed, and the trend 
from a linear regression shown as the red line. EIE, which differs from traditional ice area, is defined as the 
total surface area, including land, north of the z onal-mean i ce edg e latitude, and thus is proportional to the 
sine of the ice edge latitude. EIE was defined bv lEisenmanI ( 1201 Oh to deal with the geometric muting of ice 
area associated with the seasonal bias of the influence of the Arctic basin land mass boundaries. Here we use 
EIE and hence refer to it either via this acronym or simply ice extent. 



2. Multifractal Temporally Weighted Detrended Fluctuation Analysis 

(MF-TWDFA) 

Here, we harness a new variant of MF-DFA to study the long-range correlations and fractal 
scaling properties of time series data of basin wide Arctic sea ice extent and satellite 
albedo retrievals. The new variant, called multifractal tempor ally weighted de t rended 



fluctuation analysis (MF-TWDFA), was recently developed by IZhou and Leung! (l2010h 



and uses weighted regression to express the intuition that points closer in time are 
more likely to be similar than those that are more distant. To insure that this paper 
is reasonably self contained we motivate the approach by first discussing the general 
rationale for multifractals versus the treatment of time series using autocoiTclations and 
spectral analysis. We then describe MF-DFA followed by a development of the aspects 
that distinguish MF-TWDFA from MF-DFA. 

(a) Rationale for a multifractal approach 

Consider a time series of length i = l, ...,N describing a quantity Xj. The linear two- 
point autocorrelation function C{s) for realizations of Xi separated by an increment s is 
given by 

N-s 

^(^) = ^^w^) E (^^ - (^^+^ - ^0 ' (2.1) 

in which a\ is the variance and Xi is the mean. Long-term persistence in a time series 



occurs when on average the linear correlations become arbitrarily long (e.g., see iFederl 

N 



19881) . Whence, C(s) oc s and the mean correlation time 7^ = Jg C{s)ds diverges 
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as ^ oo for < 7 < 1- Whereas, when 7 > 1 we have finite Ts and hence we say 
that the Xi are short-term correlated. Finally, when C{s) = for s > then the Xi are 
uncorrelated. Therefore, there may be a time s* that delimits correlated from uncorrected 
Xi viz., C(s < s*) > and C{s > s*) = 0. 

Despite the simplicity of characterizing a system using a single scaling exponent 7 
it has long been understood that a wide range of natural a nd laboratory system s have 



correlations that cannot be captured by such an approach dStanley and Meakinl . Il988 



Fedeil.[l988h . It is particularly important to extract the scaling behavior for the long-term 



correlations, wherein the variability on long time scales may be poorly represented by a 
single scaling exponent and long-term trends or periodicity may obscure the calculation 
of the autocorrelation function. Moreover, the distinction between trends and long-term 
correlations in stationary time series can become blurred. These comprise some of the 
reasons for the introduction of a multifractal description of the correlations which replaces 
a single value with a continuous spectrum of scaling exponents. 

(b) Multifractal detrended fluctuation analysis 
There are four stages in the implementation of MF-DFA dKantelhardt et all l2002h . 



Firstly, one constructs a nonstationary profile Y{i) of the original time series Xi, which 
is the cumulative sum 

i 

Y{^ = Y,{Xk-Xk), where i = l,...,N. (2.2) 

fc=i 

Secondly, the profile is divided into Ns = int(A^/s) segments of equal length s that do not 
overlap. Excepting rare circumstances, the original time series is not an exact multiple of 
s leaving excess segments of Y{i). These are dealt with by repeating the procedure from 
the end of the profile and returning to the beginning and hence creating 2Ns segments. 
Thirdly, within each of the z/ = 1, ...2Ns segments the variance Vai-(z^, s) of the profile 
relative to a local least squares polynomial fit yv{i) of order is 

1 " 

Var(i/, s) = -y] [Y{[u - l]s + i) - y,{i)f . (2.3) 
Finally, the generalized fluctuation function is formed as 



(2.4) 



and the principal tool of MF-DFAn is to examine how Fq{s) depends on the choice of 
time segment s for a given degree of polynomial fit n and the order q of the moment taken. 
The scaling of the generalized fluctuation function is characterized by a generalized Hurst 
exponent h{q) viz., 

(2.5) 

When the time series is monofractal then h{q) is independent of q and is thus equivalent 
to the classical Hurst exp onent H. For the case of g = 2, MF-DFA and DFA are equivalent 



dKantelhardt et al.l l2002h . Whence, a time series with long-term persistence has h{2) 
1 — 7/2 for < 7 < 1. However, short-term correlated data, decaying faster than 1/s, has 
7 > 1 and finite Ts leading to a change at s = s*, and asymptotic behavior defined by 
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h{2) = 1/2. Moreover, the connection between h{2) and the decay of the power spectrum 
S{f) oc f^^, with frequency / is h{2) = (1 + /3)/2 (e.g., iRangaraian and Dind . l2000h . 



Therefore, one sees that for classical white noise, /3 = and hence h{2) = 1/2, whereas 
for Brownian or red noise /3 = 2 and h{2) = 3/2. 

In general multifractal time series exhibit a scaling dependence on the moment q. One 
can relate the g eneralized Hu rst exponents that we focus on here to other multifractal 



exponents (e.g., iFederl 119881) and this is the subject of a sepai^ate publication. Such 
complimentary exponents are useful because (a) different exponents are more easily 
extracted from different data sets and (b) multifractality can originate both in a broad 
probability density as well as large and small fluctuations having a different long-terrn 
persis tence. The MF-DFA procedure can distinguish between these ('Ka ntelhardt et al. 



l2002h and thus different exponents can provide tests of distinct multifractal origins. 
Here, in Section |^ we focus on the generalized Hurst exponents, which make the most 
transparent contact with the observations that are the topic of our study. 

(c) Moving Windows and MF-TWDFA 

The generalization of MF-DFA by IZhou and LeungI (|2010|) applies a variant of the 
weighted least squares approach to fitting the polynomial y,^{i) to the profile Y{i) on 
each interval u. Here, rather than using nth order yi/(i)'s to estimate Y{i) within a fixed 
window, without reference to points in the profile outside that window, a moving window 
which is smaller than s but determined by distance between points is used to construct a 
point by point approximation yu{i) to the profile. Thus, instead of Eq. l2.3l we compute the 
variance up (u = 1, ...jNg) and down (u = Ng + 1, 2A'^s) the profile as 

Var(z., s) = i Ei=i{>^([^ - 1]^ + vii^ " 1]^ + i)V 
for v = 1, Ng, and 

(2.6) 

Var(z., s) = i ELiiyi^ -W- Ns]s + i) - y{N -\v- iV,]s + i)f 

forz^ = iVs + l,---,2iVs. 

Therefore we replace the global linear regression of fitting the polynomial yv(i) to the data 
as appears in Eq. (12. 3I ). with a weighted local estimate yvii) determined by the proximity 
of points i to the point i in the time series such that |i — j| < s. A larger (or smaller) 
weight Wij is given to ?/j^(i) according to whether |z — j| is small (large). 

In the ordinary least squares method one minimizes the sum of the squared differences 
between the predicted and the original function, whereas in weighted least squares a 
weight factor is applied to the squared differences before performing the minimization. 
There are many varieties of weighted l east s quares originating in time series analysis for 



evenly spaced observa tions (iMacaulevL 



1931h and "robust" locally weighted regression for 



unevenly spaced data (IClevelandl. 1 1979 ). The standard weighted linear regression scheme 



fits observed variables Yii) and X^^ii) assuming a Unear relationship as 

Y(i) = ao{i) + ^ak{i)Xk{i) + (2-V) 



k=l 
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where e{i) is the error, which may (but need not) be homoscedastic. Here, ak{i) is the fc* 
fitting parameter at time i, the vector a of which is determined from 

a={X'^wX)-^X'^wY, (2.8) 

where w is a diagonal matrix with weight elements Wij with magnitudes that depend 
on the proximity to i thereby estimating an afc(f) according to a prescribed criterion of 
proximity. Hence, y,y{i) is determined by a weighted least squares approach using 

y^{i)=aoii)+ai{i)i + e{i), i = l,...,N. (2.9) 

Here again the fitting pai^ameters a{i) = [ao(i), depend on the place i in the 



series and are determined from equation (12.81) in which the Wij, with j = 1, are 
the diagonals of an n x n diagonal matrix, and A' is an n x 2 matrix with the first column 
being unity and the second running from 1 to n. We use the same bisquare function for 
the weights 









s 



2X 2 



wu = i \ L ^- J y ' , (2.10) 

otherwise, 

as did IZhou and Leung! (l2010h . although one can envisage a range of possibilities for 



dealing with the temporal weighting. 

We conclude this section by emphasizing that a principal advantage of MF-TWDFA 
over conventional MF-DFA is the robustness of extracting the scaling of the fluctuation 
function, and hence the crossover points between scalings, that indicate the underlying 
processes reflected in the data. Of particular importance is the fidelity of extracting the 
long-term scaling behavior in geophysical data sets which can often be obscured in MF- 
DFA and, as we shall see below, compromises MF-TWDFA when a strong seasonal cycle 
remains in the original time series. Moreover, in MF-DFA the profile of the time series 
is fit using discontinuous polynomials, which can introduce errors in the determination of 
crossover times for new scalings, and can be particularly questionable at long time scales. 
Finally, for time series of length A^, whilst MF-DFA is typically informative only up to 
N/A, MF-TWDFA can be canied out to N/2. 



3. Sea Ice Geophysical Data 

We use MF-TWDFA to examine the multi-scale structure of two satellite based 
geophysical data sets for Arctic sea ice; the Equivalent Ice Extent (EIE) and albedo 
retrievals from the Advanced Very High Resolution Radiometer (AVHRR) Polar 
Pathfinder (APP) archive. The EIE data derives from retrievals of satellite passive 
microwave radiances over the Arctic converted to daily sea ice conc entration u s ing th e 
NASA Team Sea Ice Algorithm in the same manner as described by lEisenma nl (I201(]h . 



who focused on the origin of the difference between ice extent and EIE. We refer the 
reader to Eisenman's paper for a detailed description of the determination of EIE. The 
mean EIE seasonal cycle from 1978-present is shown in Fig.[2l Daily satellite retrievals 
of the directional - hemispheric apparent albedo are determined from the APP archive as 
descr ibed in a separate pub lication dedicated to a different form of analysis than presented 
here dAgarwal et al.L 1201 ih . The apparent albedo is what would be measured by upward 



and downward looking radiometers and thus varies with the state of the atmosphere and 
the solar zenith angle. 
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Figure 2. The mean seasonal cycle of the Equivalent Ice Extent (EIE) during the satellite era (Fig.[T]( shown 
as the solid line, with the dashed lines denoting one standard deviation. 
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Figure 3. Albedo histograms shown for days in mid-March (black) and -September (white). If there is ice in 
a pixel for the 23 year record then we compute the albedo for that pixel. We com pute a histogram for each 
day of the year from the number of pixels for each average albedo bin for that day jAgarwal et alll20lll) . 



The APP dataset ha s been refined for us e in a wide range of polar studies and 

and refs. therein). In brief, the AVHRR 



is described in detail in lAgarwal et al.l (12011 



channels range from the visible to the thermal infrared (0.58 - 12.5 fim) and measure top 
of the atmosphere reflectances and brightness temperatures. With 5 km x 5 km resolution, 
we analyze albedo retrievals from 1 January 1982 through 31 December 2004, taken 
daily at 1400 hours. Sea ice is distinguished from land and open water using microwave 
brightness temperatures and filtering the Surface Type Mask data with the NASA Team 
Sea Ice Algorithm which distinguishes between first-year (FYI) and multi-year ice (MYI) 
concentrations. The approach ascribes a MYI flag to a region containing at least 50 percent 
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Figure 4. (a) The profile for the EIE and (b) albedo using Eg. 12. 21 with a seasonal cycle in the original time 
series, (c) The profile for the EIE and (d) albedo using Ea. \2.2\ removins the seasonal cycle from the original 
time series. 

of this ice type, with uncertainties depending on the season (e.g., melt pond fraction) and 
region (near the ice edge), along with the surface type categories. On physical grounds 
the albedo data are filtered to remove any values greater that 1 or less than 0.2. Each pixel 
is assess ed eve r y day for the presence of ice and then albedo is averaged for that pixel 
( Agarwa l et al.l llOllI) . Examples of histograms for mid-March and mid-September are 
shown in Fig. [3j 



4. Results and Discussion 

Here we focus most of our discussion on the nature of the fluctuation functions Fq{s) that 
emerge from MF-TWDFA and return later to the exponents, generalized dimensions and 
singularit y spectra. First, we note that we have eliminated the possibility of spurious multi- 
fractility dSchumann and K antelhardt. 2011) by checking the usual measures A/i2o = 
20) — /i(20) and Aa. Second, whilst we show the fluctuation functions for a range 
of q, for clarity we also show plots solely for q = 2 and remind the reader that for 
h{2) = 1/2 the system exhibits a completely uncorrelated white noise dynamics, and for 
h{2) = 3/2 there ai^e red (or Brownian) noise coiTclations. Whereas for < h{2) < 1/2 
the dynamics are anticoiTclated. Moreover, as noted above, for analysis by power spectra 
/3 = 2/i(2) — 1. It is evident from equation (12.41) that small temporal fluctuations are 
characterized by h{q)'s for g < whereas large temporal fluctuations are characterized 
by /i(g)'s for q>0. As we discuss below the fluctuation functions Fq{s) do indeed 
distinguish the scaling of small and large fluctuations. 
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Figure 5. (a). The fluctuation function from equation i2.5\ for the Equivalent Ice Extent (EIE) with a seasonal 
cycle. The q's are shown in the panel and s is measured in days throughout, (b). The fluctuation function from 
equation J2.5b for Equivalent Ice Extent (EIE) with a seasonal cycle for q=2. The stars denote the crossover 
times associated with a slope change at approximately 13 days, 1 year and 5.25 years. The blue line and the 
red segment of the black curve both denote white noise with h{2) = 1/2. (c) and (d) correspond to (a) and (b) 
but without a seasonal cycle. Here, the stars denote crossover times of approximately 8 days, 1 year, 2 years, 
6.9 years and 8.9 years. Again, the blue line and the red segment of the black curve both denote white noise 
with h{2) = 1/2. 
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Figure 6. (a). The fluctuation function from equation \2.5l for tfie ice albedo with a seasonal cycle. The q's 
are shown in the panel, (b). The fluctuation function from equation l l2.5t for ice albedo with a seasonal cycle 
for q=2. The stars denote the crossover times associated with a slope change at approximately 7 days, 16 
days, 1 year and 2.2 years. The blue line and the red segment of the black curve both denote white noise with 
h{2) = 1/2. (c) and (d) correspond to (a) and (b) but without a seasonal cycle. The stars denote crossover 
times of approximately 8 days, 17 days, 487 days, 7.5 years and 9 years. Here again, the blue line and the red 
segment of the black curve both denote white noise with h{2) = 1/2. 
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On all segments of the Fq{s) plots the slopes are computed and the "crossovers" 
refer to the times s where the curve changes slope. As a test of the fidelity of the many 
crossovers in slope (associated with particular instrinsic time scales) detected with MF- 
TWDFA, we applied both MF-DFA and MF-TWDFA. We found that the crossovers for 
time scales of 2 years or longer would not have been captured by MF-DFA because of 
large amplitude fluctuations in this range. Note, for the EIE data with a seasonal cycle this 
result was confirmed for polynomial fitting in the detrending step of MF-DFA of up to 
order 9, while for the other data up to order 3. Moreover, and that which we focus on here, 
the MF-TWDFA analysis used on the time series after removing the seasonal cycle leads 
to the extraction of crossovers associated with long term persistence that are "masked" 
when the seasonal cycle is not removed. Next we discuss this finding in more detail. 

(a) Masking long term correlations by strong seasonal cycle 

A prominent feature of our analysis is the role played by the seasonal cycle. We lay 
bare the distinction between profiles with and without the seasonal cycle in Figs. |4l which 
is the origin of the striking distinction between the upper and lower panels of Figs. [5] 
and [6] for the EIE and albedo respectively. It is seen that the strength of the seasonal 
cycle is such that it "masks" the dynamics on time scales longer than seasonal thereby 
suppressing the fingerprints of long term persistence. This is clearly displayed in the 
fluctuation functions. Firstly, in Figs. |5] and [6] panels (a) and (b) which have a seasonal 
cycle, we see that the curves for all g's converge at the 1 year time scale. This convergence 
is removed when the seasonal cycle is removed as seen in Figs. [5] and [6] panels (c) and 
(d). Secondly, when the original time series still possesses the seasonal cycle, the slope 
abruptly drops below h{2) - 1/2 at the 1 year time scale to transition to an anticorrelated 
structure. This behavior is distinct from the clear transitions and positive slopes at longer 
time scales that are seen when the seasonal cycle is removed from the original time series. 
Whilst there is a finite union of underlying processes that influenc e the ice exten t and the 
ice albedo, that union is cleai^ly not one to one (see discussion in lAgarwal et al... .201 1.) . 
Thus, although we find similar qualitative structure in the analysis of the EIE and albedo 
data with and without a seasonal cycle, there are important quantitative distinctions. For 
example, the high frequency weather time scales of order week's are clearly seen in all of 
the data. However, Figs. |5] and [6] (a) and (b) clearly show masking beginning at one year, 
the transition to anticorrelated behavior and then a steeper-white noise-slope of h{2) = 
1/2 returning at 5.25 and 2.2 years in the EIE and albedo respectively. Whence, all long 
term persistence is destroyed by the strength of the seasonal cycle which is associated 
with the strong periodicity remaining in the profiles of Fig. |4l We see that the general 
steep/flat/steep behavior of the slopes of the fluctuation function with increasing time is 
seen in all panels of Figs. [5] and [6l However, the essential distinction is that when the 
seasonal cycle is removed long term persistence is reentrant. Indeed, time scales of 7 
and 9 years aie revealed in both the EIE and albedo data. 

{b) White noise on seasonal time scales 

The region between 1 and 2 years in the EIE data without a seasonal cycle shows 
h{2) 1/2. In other words, although the strength of the seasonal cycle is such that it 
dominates the power spectrum (not shown), when removing the seasonal scale from the 
original EIE data the system exhibits a white noise behavior from the seasonal to the bi- 
seasonal time scales. However, the clear fingerprints of the short (weather) and long (~ 
7 and 9 year) time scales remain. The implications of this finding are, unfortunately for 
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the goal of forecasting, rather clear; a given ice minimum (maximum) can be followed by 
a minimum (maximum) with a larger or smaller magnitude. The same logic follows for 
minima/maxima separated by two years and the EIE states on time sea les between. Most 
of the discussion in the literature focuses on the extremes in the observations (maxima 
and minima) which we find here to exhibit a nearly white time distribution. If we interpret 
the results in the limit h{2) = 1/2 then the autocoiTclation will have strong peaks at the 
seasonal maxima, minima and bi-seasonal maxima and minima with rapid decays for 
times from one to two years. As discussed below, this is consistent with aspects of other 
studies in which rapid decoiTclations ai^e found. 

It is important to be clear regarding what these results do not tell us. We are not saying 
that existence of minima or maxima (or states between) are uncorrected in time, but 
solely that the magnitudes of those states cannot be predicted to be larger or smaller 
from these data alone. Moreover, the lack of an extended autocorrelation does not a-priori 
constrain the probability distribution(s) associated with the process(es) from which the 
value of the observable originates, but rather solely refers to the temporal distribution of 
the observable. In other words whilst a particular- probability distribution may underly 
the nature of how a given magnitude of the EIE or albedo is reached, we find that the 
temporal distribution of the sequence of observed magnitudes is white on seasonal to bi- 
seasonal time scales. Thus, on such time scales these data are short-term correlated, are 
associated with a finite Ts and an autocoiTclation that decays faster than ~ s~^. Therefore, 
on these short time scales, the processes controlling the EIE or albedo states in the Arctic 
are the familiar regular strong seasonal scale radiative and advective forcings. Indeed, this 
is explained by the robustness of the response times associated with the seasonal control 
of the ice state through the competiti on between the ice-albedo feedback in summer an d 
the long-wave forcing in winter (e.g.. iMoon and Wettlaufed . 1201 ll : iTietsche et all 1201 ih . 
However, on longer time scales, we can associate the retreat of the ice cover with the 
observation of the upswing in h{2) towards a more "Brownian" dynamics; h{2) ~3/2. 
Such trends, when accumulating future data over decadal time scales could well be 
compromised by a higher degree of noise in the state of the system, but to speculate 
more moves beyond the level of firm evidence. 

We test the idea that the rapid decline in the ice coverage during the most recent decade 
is dominating the longest crossovers found in our analysis as follows. The EIE record 
is approximately a decade longer than that for the albedo and covers the decade of the 
most rapid decline. Thus we ask what crossovers are seen when we reanalyze this record 
between 1982-2004, the time range corresponding to the albedo data set? It is found that 
the two long term crossovers decrease from 7 to 6.2 and from 9 to 8.4 years. Moreover, 
the seasonal to bi-seasonality of the white noise persists. Therefore, the increase in the 
crossover times as we include the recent data appeai^s to be the most plausible reflection 
of the recent declines. However, we leave for future study what would likely be a less 
quantitative endeavor of conelating these longer time scales with climate indices that 
have time scales ranging from 2 to 7 years (for the El Nino Oscillation) to a decade or 
many decades (such as for the Pacific Decadal, North Atlantic or Arctic Oscillations). 



(c) Generalized Hurst exponents 

Figure |7] shows the trends in the generalized Hurst exponents, h{q), associated with 
the strong seasonal cycle. We calculate these exponents, which describe the slope of 
the fluctuation curve, for each q as the linear approximation of the entire curve. It 
is principally used here to demonstrate the strong influence of the seasonal cycle in 
suppressing h{q) for all q. There are several important effects immediately evident. First, 
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Figure 7. The generalized Hurst exponents as functions of q for Equivalent Ice Extent (EIE) (a) with a 
seasonal cycle and (b) without a seasonal cycle, and the ice albedo (c) with a seasonal cycle and (d) without 
a seasonal cycle. 
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the general reduction in the magnitudes of the h{q) associated with the seasonal cycle is 
synonymous with the suppression of long term persistence, restricting the crossovers to 
short term dynamics and masking behavior on time scales longer than seasonal, as was 
seen in the fluctuation functions themselves in Figs. |5] and [6l Second, this suppression 
of long term persistence by definition suppresses multiple crossovers as well as trends on 
long time scales. Third, the combination of the reduced magnitude and the general flatness 
of the curves with the seasonal cycle present indicates a suppression of the "dynamic 
range" of multifractality. The analysis and interpretation of related multifractal exponents 
is the subject of a separate publication. The main point of Fig.|7]is to display the behavior 
and interpretation of Figs. [5] and [6] in a more compact manner. 



(d) Discussion 

In his EIE observational analysis lEisenman (I201(]h noted that the overall retreat 
of the ice cover is governed by a noisy signal and particulai^ly so at the September 
minimum. Observations using 30 year s of monthly ice extent and area data show 
rapid decorrelations of a few months (iBlanchard-Wrigglesworth et all. 1201 ll) using lag 
correlation method assuming that the underlying process is a first-order autoregressive 
process. A similar approach was applied to fit hindcasting si mulation output whic h 
showed a persistence time scale for September ice area of 1 2 vears ^Armour et al.Ll201 lah . 
In a study of the forecast skill of a linear empirical model iLindsay et all (l2008h found no 
skill in predicting detrended data for time scales of three or more mont hs. Finally, the 



rapid d ecay in ice area anomalies is also found in climate models (see lArmour et al 



201 Ibl. and refs therein) and this is argued to underlie the lack of hysteresis associated 



w ith the loss of su mmer sea ice in both global models and theoretical treatments such as 



in 



Eisenman and Wettla ufer ( 2009.) 



The corpus of these studies lead to the general conclusion of rapidly decaying 
correlations and hence compromised predictability. Thus, in this respect they are 
heuristically consistent with our interpretation here. However, there remain important 
distinctions. Firstly, even a single moment of the fluctuation function demonstrates the 
existence of multiple time scales in the data which, as noted in Section |^ cannot 
be treated in a quantitatively consistent manner with a single decay autocoiTclation. 
Moreover, we showed in Section |^ that if the seasonal cycle is not removed one will 
always observe a single crossover time (longer than the synoptic time scale) of 1 year; 
a time scale at which all moments converge. Thus, the upper bound on the persistence 
time in any study that assu mes an autoregressive process will inev i tably be ~ 1 year, as i s 
indeed found for ice area dBlanchard-Wrigglesworth et al.l . 1201 ll : lArmour et al.l 1201 lah . 
Apart from the single moment evidence of multiple scales in the system, the approach 
here highlights the dangers of not carefully detrending the seasonality and dealing with 
stationarity. Secondly, both data sets as used for this analysis and many of the studies 
mentioned above do not explicitly distinguish between ice types. Thus, this analysis 
applies to the aggregate of the ice cover, whereas, subject to a wide range of caveats, 
many models can diagnose features such as thickness or type. Present work focuses on 
extending this analysis in to various ice types. 



5. Concluding Remarks 



We have examined the long-term correlations and multifractal properties of daily satellite 
retrievals of Arctic sea ice albedo and equivalent ice extent (EIE), for periods of ~ 23 
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years and 32 years respectively, with and without the seasonal cycle removed. A recent 
development called Multifractal Temporally Weighted Detrended Fluctuation Analysis 
(MF-TWDFA), which exploits the intuition that in any time series points closer in time 
ai^e more likely to be related than distant points, was adapted for use with these data. 
Points in the records nearer each other are weighted more than those farther away in 
order to determine a polynomial used to fit the time series profile, which is a cumulative 
sum of the time series. As a methodology the approach offers several advantages over 
the more generally applied MF-DFA. The profile in MF-DFA is fit using discontinuous 
polynomials, which can introduce eiTors in the determination of crossover times for new 
scalings, and can be particularly questionable for long time scales. Whereas MF-DFA 
is typically informative only up to N/A for time series of length N, MF-TWDFA can 
be carried out to N/2. Additionally, the generalized fluctuation functions Fq{s) for all 
moments g as a function of time scale s are substantially smoother for all s and this is 
particularly so for large values. This facilitates clear extraction of crossover times from 
one scale to another. 

The generalized Hurst exponents and multiple crossover timescales were found to 
range from the synoptic or weather time scale to decadal, with several between. Such 
multiple time scales were exhibited in both data sets and hence the approach provides a 
framework to examine ice dynamical and thermodynamical responses to climate forcing 
that goes beyond treatments that assume a process involving a single autocoiTclation 
decay, such as a first-order autoregressive process. Indeed, the method shows that single 
decay autocorrelations cannot be meaningfully fitted to these geophysical observations. 
One of our most important findings is that the strength of the seasonal cycle is such that it 
dominates the power spectrum and "masks" long term correlations on time scales beyond 
seasonal. When detrending the seasonality from the original record, the EIE data exhibits 
a white noise behavior from seasonal to bi-seasonal time scales, but the clear fingerprints 
of the short (weather) and long (~ 7 and 9 year) time scales remain, demonstrating a 
reentrant long-term persistence. Therefore, it is not possible to distinguish whether a 
given EIE minimum (maximum) will be followed by a minimum (maximum) that is larger 
or smaller. This means that while it is tempting to use an anomalous excursion associated 
with a low ice year to predict the following year's minimum, or that two year's henceforth, 
the present data do not justify such a prediction. 

We tested the idea that the rapid decline in the ice coverage during the most recent 
decade is associated with the longest crossovers found in our analysis of the EIE. This was 
done by reanalyzing this record between 1982-2004, the time range corresponding to the 
albedo data set. In so doing, we found that the exclusion of the most recent decade led to 
the two long term crossovers decreasing from 7 to 6.2 and from 9 to 8.4 years. Moreover, 
the white noise structure on seasonal to bi-seasonal times is the same in both records. 
Hence, the increase in the crossover times associated with including the last decade of 
data appears to be the most plausible reflection of the recent declines. However, we have 
not attempted to correlate these longer time scales with other climate indices. 

Finally, other methods of analyzing such data find solely a rapid decorrelation, whereas 
we find multi-year and decadal transitions as well as the origin of the dominance of 
the seasonal cycle in long term persistence. Hence, we believe that combining such 
multifractal studies of model output and other observations will substantially improve 
the acuity with which one can disentangle the strength of the seasonal cycle in this highly 
forced system from the longer term trends. 
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